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Solubility isotherms for CaHP0 4 -2H 2 0, dicalcium phosphate dihydrate, DCPD, in the 
ternary system Ca(OH) 2 — H. ? PO»— H 2 were determined at 5, 15, 25, and 37.5 °C in the pH range 
3.5—7; the relative positions of the isotherms indicate that DCPD has a negative thermal coefficient 
of solubility. The solubility product, /C,, of DCPD and the stability constants K x and K u for the ion pairs 
[CaHPO^] and [CaHjPOj ], respectively, were obtained as functions of temperature by the use of a 
generalized least squares procedure subject to three condition functions — constancy of the solubility 
product, electrical neutrality in the solution, and congruent dissolution of the solid. The equations 
obtained are 

In K s = - 8403.5/7+ 41 .863 - 0.096787' 
In /C. r =51090/7 , -341.14 + 0.5880r 
In K 9 = 19373/7 1 - 122.81+0.19947. 

The existence of a maximum in K s in the neighborhood of 25 °C is plausible on the basis of available 
thermodynamic data for DCPD. Thermodynamic functions are reported for the solution of DCPD and 
for the association of the ion pairs. 

Key words: Calcium phosphate; dicalcium phosphate dihydrate; ion pairs; solubility; solubility 
isotherms; solubility product. 



1. Introduction 

In a previous publication [l] 1 on the solubility at 
37.5 °C of dicalcium phosphate dihydrate (DCPD), 
CaHPC>4 * 2H 2 0, in dilute phosphoric acid solutions, 
we reported that the "apparent" solubility product, 
K f DC , exhibited a concentration dependence. This de- 
pendence could be removed by postulating the pres- 
ence of the ion pairs [CaHPOSJ and [CaH 2 POj], 
making it possible to calculate a "true" solubility 
product, as well as the stability constants of these ion 
pairs. 

In view of the importance of DCPD in the areas of 
biology, agronomy and industry, we have extended to 
"other temperatures our solubility measurements of 
this compound in the ternary system Ca(OH) 2 — H 3 P0 4 
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— H 2 to obtain standard free energies, enthalpies and 
entropies for dissolution of DCPD and for formation 
of the ion pairs. 

We report here the consolidated results and calcula- 
tions on the solubility behavior of DCPD at 5, 15, 25, 
and 37.5 °C. Although measurements at 25 °C had 
been reported previously [2], it was decided to repeat 
the work at this temperature in order to make all the 
data comparable in terms of analytical techniques. In 
the light of our continued experience with these sys- 
tems we have reevaluated the data collected at 37.5 
°C [1] by the use of a more general statistical treat- 
ment. No measurements formerly reported [1] have 
been rejected, but additional measurements have been 
included. 

The data were treated by generalized least squares 
[3, 4] to obtain estimates of the constants of interest, 
their standard errors and the adjusted values of the 
data. The major departure from customary [5, 6] least 
squares applications introduced in this paper is the use 
of multicomponent observation and condition function 
vectors: each of TV independent systems is associated 
with a 4-component observation vector and a 3- 
component function vector. It will be shown that, 
where applicable, the imposition of more than one 
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condition function per system materially reduces the 
computed standard deviation. We have applied multi- 
function conditions in other solubility studies [7] and 
in the study of diffusion in ternary systems [8], 



2. Methods and Materials 



2.1. Preparation of Materials 



Dicalcium phosphate dihydrate and the phosphoric 
acid solutions were prepared and characterized as 
described in previous publications [1, 2]. 

2.2. Experimental Method 

The DCPD preparations were leached in a thermo- 
statted saturator with a series of dilute phosphoric 
acid solutions as described previously [1]. The pH of 
the effluent was measured in situ and a number of 
samples (4—8) were collected for the determination of 
total calcium and phosphorus. 

2.3. Calcium and Phosphorus Determinations 

Phosphorus was determined spectrophotometrically 
(relative error, 0.01) in a double beam instrument using 
the vanado-molybdate reagent of Brabson et al. [9], as 
described for our previous work at 37.5 °C [1]. Calcium 
was determined by atomic absorption with a relative 
error of 0.015. 

2.4. pH Determinations 

The pH measurements were made with glass and 
calomel electrodes with a commercially available pH 
meter (claimed relative accuracy ± 0.0037 pH unit). 
Standardization of the meter was done with certified 
NBS buffer standards [10]. The error in the pH meas- 
urements was taken as 0.015, which includes the esti- 
mated error (0.005) in the value of the standard buffer 
used [10]. 

3. Calculations 

The aim of these calculations was to obtain as a 
function of temperature, the "true" solubility product, 
K s , of DCPD, the stability constants K. x and K y , of the 
ion pairs [CaHPO*}] and [CaP^POj ], respectively, and 
the adjusted values of the experimental data. This was 
accomplished as described, by minimizing the weighted 
sum of squares of the residuals in the observations 
subject to three condition functions: saturation with 
respect to DCPD, electroneutrality, and congruent 
dissolution of this salt. 

Preliminary calculations were done separately on 
the data collected at each of the four temperatures and 
the equilibrium constants, K m , so obtained were each 
fitted by subsequent adjustments to equations of the 
form In Km = A mi T~ 1 + A m2 + A m . i T. Since the results 
of this procedure showed that the assumed form of 



temperature dependence was quite reasonable, the 
calculations reported here were then done on the 
combined sets of raw data using the above form of 
temperature dependence ab initio. In this way the 
nine constants were obtained without the necessity 
of further adjustment. Throughout this work the 
temperature was considered as a preset, error-free 
variable. 

In order to formulate the condition functions, it 
is first necessary to characterize the systems in terms 
of the measurements and constants. This is done in 
sections 3.1, 3.2, and 3.3. 

3.1. Balance Equations 

Using the symbols X° and Y + for the names and con- 
centrations of the species [CarlPOJJ] and [CaH 2 POj], 
respectively, the total molar concentrations of calcium, 
Ca, and phosphorus, P, are given by 



Ca=[Ca ++ ]+r 



P = J[H,POi- 3 ]+T 



(1) 



(2) 



in which t = X°+Y + , and [Ca ++ ] denotes the concen- 
tration of free calcium ions. Equation (2) may be written 
in terms of the activity of the HPO4 ion, (HPO4 ), and 
the three dissociation constants k\ , k 2 , and k 3 of 
phosphoric acid as 



where 



N= J N igi 



P=(HPOj) TV-hr, 



(H+) 2 (H+) £ 3 



(3) 



kik 2 



(H+) 



#3. 



(4) 



In eq (4), (H + )=10 _pH is hydrogen ion activity and 
gi is the reciprocal of the activity coefficient fi of 
species i. The subscript i takes the value of the charge 
magnitude of the several free phosphate species. 

3.2. Solubility Product and Concentrations of Ion 

Pairs 

The solubility product is defined here as 

K D c=(Ca++)(HPOr)=(Ca-T)(P-T)/«r Ca iV. (5) 

In eq (5), and in eqs (6, 7) below, the parentheses 
around quantities with a charge symbol denote activi- 
ties. The third member of eq (5) follows from eqs (1) 
and (3). An apparent solubility product, K^c, was 
calculated taking r to be zero. K DC , a $ defined by eq (5), 
is distinct from the quantity K s ; the latter is the least 
squares value of the solubility product constant as 
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found by the adjustment procedure described in 
section 3.5. The quantity K DC , on the other hand, is 
the ionic-activity product computed for each indi- 
vidual system; in the absence of experimental errors, 
its value should be the same for all systems, and equal 
to K s , if the saturation condition is fulfilled and if the 
model for ion pairs is valid. 

The concentrations of the ion pairs can be expressed 
in terms of the stability constants as follows: 

K x = (CaHP02)/(Ca++) (HPO=) =Z /£ DC , (6) 

K, = (CaH 2 PO + )/(Ca++)(H 2 P0 4 -) = Y^K^giNi. (7) 

Adding eqs (6) and (7) and using the identity t=X°-\-Y + , 
a quadratic equation in r is obtained whose pertinent 
solution is 



r=U-(U 2 -P-Ca) 1 



/2 



(8) 



where E/ = i{P + 0* + g c *N I {K x + gtNiKy)}. Writing 
R = (1 -\- Kxlg\NiKy)~ l , the ion pair concentrations 
are 

A° = (I-R)t, Y+ = Rt. (9) 

3.3. Ionic Strength and Activity Coefficients 

The ionic strength, /x, is defined in these calculations 
by 

2 M = 4(Ca-T)+£ H (H + )+£oH^/(H + ) + F + 

+ <MP-t)//V, (10) 
where k w is the dissociation constant of water, 

i = 

and the JVi are defined by eq (4). 

Inasmuch as individual ionic activities are required 
for the computations, the results are bound to reflect 
the particular expression arbitrarily chosen for the 
activity coefficients. For this reason calculations were 
carried through using two different expressions: the 
Debye-Huckel limiting law extension (D.-H.) and the 
expression used by Davies [11]. These expressions 
are given by eqs (11a) and (lib), respectively. 

-log/i=/i(r)zfVM/(l-hS(r)aiVii) (11a) 

-log fi=A(T)zf{ Vm/0 + V£) -0.3**} ( llb ) 

where z\ is the valence of the ah ion, a, is the "distance 
of closest approach" parameter, and A and B are 
temperature-dependent constants. Following conven- 
tion, the activity coefficients of the uncharged species 
H3PO4 and [CaHPO<j] were taken as unity. The 
coefficient of the charged ion pair [CaH^POl] was 
taken to be equal to that of the ion H 2 POj. 



Values used at the various temperatures for the 
dissociation constants k u &2, & 3 , kw, and for the 
parameters A(T) and B(T) are given in table Al, 
section 6.2. In table A2 are listed the values for the 
parameters a\. 

Iteration Procedure. Approximate values were first 
obtained for N in eq (4), and for r and F + in eqs (8) 
and (9), using initial guesses for the constants K x and 
Ky and assuming all//=l. An approximate value of 
the ionic strength, /x, eq (10), was then obtained which 
yielded better estimates for the/} (eq (11)). The calcula- 
tions were iterated until a constant value of jjl was 
attained (2-3 cycles). New estimates for the parameters 
A m k were then obtained by the adjustment procedure 
to be described in 3.5; new values of K x and K y com- 
puted from eq (16), when inserted into eqs (8) and 
(9), yielded improved estimates for the ion pair 
concentrations. A new iteration on /a was started, but 
using the last values of the//, and the whole process 
was continued until the weighted sum of squared 
residuals attained a constant value. 

3.4. Weights of the Observations 

The weights of the observations, defined below by 
(12), were computed from the estimated variances. The 
constant So was taken as 1.0 X 10 _H . 



w = sf) /variance 



(12) 



Using the value given in section 2.4 for the error, the 
variance of pH was calculated as (0.015) 2 . The vari- 
ances of phosphorus and of calcium were computed 
from the relative errors given in section 2.3. At 37.5 
°C, the relative error for calcium was 0.025. 

3.5. Least Squares Adjustment 

This section adheres closely to the notation of 
section 11.3 of reference [4]. Equation numbers re- 
ferred to in this text from that source are enclosed in 
square brackets. The symbols used in this paper are 
defined in table 1 ; they differ from those in table [1 1.3.3] 
of reference [4] in several respects, as indicated in 
section 6.1. 

Table 1. Matrix symbols 



Name 


Symbol 


Order 


Element 


Condition function 


fa 

yCt 

v a 

U a 

W a 

6 
J a 

fa 

L" 

A" 

c 


rXl 

nX 1 
71 X 1 
nX 1 
n X n 
pX 1 
pX 1 
rXp 
rXn 
rXr 
rX 1 
pXp 


F;(Yf; ft) 
Y? 


Observation 


Residual 


Vf 


Adjusted observation 

Weight 


Uf 

Wfj (diagonal) 


Parameter 


Correction 


Lagrangian 


F&=dF<>ldp v 
F$=dF*ldYf 

^qs ^' r qi r si w ii 
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The general formulation of reference [4], applied 
to the combined data (35 points) with 3 condition func- 
tions per system and 9 parameters, leads to the matrix 
[11.3.3, la] of order 3X35 + 9=114. In an alternative 
formulation of the normal equations [11.3.3.1, 5a], the 
matrix L is of order 3 X 35 = 105. However, if the condi- 
tion functions for the ath system depend on the meas- 
urements on that system only, as is true here, then, as 
shown in section 3.5.3, the order of the largest matrix 
that appears is 9, the maximum number of parameters 
considered here. 

a. Condition Functions 

The available conditions to be imposed on each 
system are: 

1. Ff: Saturation with respect to DCPD 

(Ca-T)(P-T)/*cJV-K,=0. (13) 

2. Ff: Electroneutrality 

2(Ca-r) +£ H (H+) +Y+-g 0H k w /(H+) 

-(/.,(P-r)/yV = 0, (14) 

3 

where </>,= V iNigi. 

3. Fg: Congruent dissolution of DCPD 

Ca/(P-P () )-1 = 0, (15) 

where P () is the molar concentration of initial phos- 
phoric acid. 

In eq (13), K s is the solubility product constant to 
be estimated. The stability constants to be estimated, 
K x and K u , appear implicitly in r and Y + ; see eqs 
(8) and (9)." 

b. Choice of Parameter and Indices 

In these calculations the equilibrium constants were 
represented by equations of the form 



lnK n 



L 

2> 



L=lor3 
01=1, . . ., c 



(16) 



where m=\, 2, 3 corresponds to the subscripts 5, x 
and y, and T is the temperature in kelvins. In the 
case L=l (computations at individual temperatures), 
T~ l was replaced by 1. The number of constants, C, 
was 3 in every case except when ion pairs were 
assumed absent. In any event, p = CL is the total 
number of parameters, Amk, to be found, and 
J3=(Ai U . . „A lL , . . .,/* C l). 

In order to establish the soundness of the data 
treatment adopted here, i.e., assumption of the 
presence of ion pairs and use of three condition 
functions, calculations were made using the following 
combinations of parameters and functions: (I) Calcula- 
tions at individual temperatures, assuming ion pairs 
present; (II) combined data at all temperatures, with 



temperature dependence for K s , but ion pairs assumed 
absent; (III) combined data, ion pairs present, and 
temperature dependence of the K m . In each of these 
cases, separate calculations were made, using as 
conditions: (A) F? only; (B) F<\ and Fj; (C) Ff , F?, and 
Fg. (When F% is omitted, the observable P is not used, 
i.e., n = 3). These combinations are tabulated below. 
The calculated quantities reported here correspond 
to Case III— C. 



Case 


c 


L 


P 


r 


n 


I 


3 


1 


3 


3 


4 


II— A 


1 


3 


3 


1 


3 


-B 








2 


4 


-C 








3 


4 


III -A 


3 


3 


9 


1 


3 


-B 








2 


4 


-C 








3 


4 



c. Reduced Normal Equations and Their Solution 

In terms of the analysis of this paper, the reduced 
or "p normal" [12] eqs [11.3.3.1, 5a], may be written 



^here 



C0 = z, (17) 



C = ^J"L«- l J", z = ^J a L«- x f«. (18) 



Since J a is of order rXp and L a is of order rXr, the 
matrix C is of order pX p. 

Equation (17) was scaled by the substitution 
= D~ 1 6* and premultiplication by D~ x to give 



C*d* 



(19) 



where D vv ,= \ 
C*=D~ 1 CD-\ 

parameter corrections is then 

0=J)- 1 C*- 



vSw' (8w' is the Kronecker delta), 
id z* = D 1 z. The solution for the 



(20) 



The inversion of C* was accomplished using Gaussian 
elimination and double precision arithmetic. 

New estimates, (3 = l3 {) — d, were thus obtained and 
used for a subsequent iteration (sec. 3.3). 

The convergence of the process at the nth iteration 
was followed by computing the sum of squares (q 2 ) 
and testing the ratio q\\q\- X against 10~ 7 . For p = 3, 
convergence was attained in 3-4 cycles and for p = 9, 
about 15 cycles were required. The expression used 
for q 1 is given in eq (A- 15) of section 6.1c. In computa- 
tions for consecutive values of a during a given cycle, 
partial sums (denoted by the prime) were accumulated: 
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C =C 4. loci a 1 la 

U (a) °(a-l) ' •/ ^ J » 

Z (a) Z (a-1)^ •/ ^ / » 



and 



« 



&=flt-n+/ b ^" 1 A 



The sum q?L=% a J**L a f a was then corrected, after 
obtaining the solution for 0, by subtracting the quantity 
SCO. 

d. Standard Deviation and Errors 

The total number of condition functions considered 
here is Nr, i.e., the number of systems, times the 
number of functions per system. Hence, [11.3.4.1, 3] 
for the estimated sum of squares is 

%(XaV a W«V«) = (Nr-p)s 2 . (21) 

The standard deviation of the adjustment is thus 

s=V[q 2 l(Nr-p)], (22) 

with dispersion (standard error) given by [13] 

<r{s}=slV[2(Nr-p)]. 

The covariance matrix of the parameter set (3 is 
s 2 C~ l by [11.3.4.2c] and (A— 1). This may be written 
as 

Cov/3 = 5 2 [(C- l )mm'],m, m' = \, . . .,c. (23) 

Here s 2 (C~ 1 ) mm >, one of the C 2 submatrices of s 2 C -1 , 
is of order L X l and is the moment matrix of the set 
{Amky Am'k')- (See eq (16)). Then the elements of the 
covariance matrix for the pair of equilibrium constants 
K m , K m > are given by 



Cov (K m , K tU ')=S 2 KmXt(C l ) mm'Xt'Kn 



(24) 



where the vector x, = (1) for L=l, and x* = 
(Ff 1 , 1, T t ) for L=3. When the selected temperatures 
are identical (t = t'), the dispersion of K m is given by 

cr{£*} = V[Cov (K m , K m )]. (25) 

For the adjusted observations the dispersions were 
computed from eq (A-10) by 



4.1. 



cr{tf?} = VTeov (tff, Uf)]. 

4. Results 

Solubility and Solubility Products 



(26) 



The compositions of saturated DCPD solutions at 
four temperatures are summarized in tables 2 through 
5. Congruent dissolution of the salt is indicated by the 



Table^. Solubility of DCPD at 5 °C 



Unadjusted quantities 





Composition of equil. solns. 




K' n( . X 10 7 


Initial 






Ca 
P-Po 






acid, P 












Mx 10 3 


pH 


Ca 
MX 10* 


P 

Mx 1() :{ 




Davies 


D.-H. 


0.0979 


7.64 2 


0.755 


0.843 


1.01 


2.96 


2.98 


.490 


6.87 8 


LOO 


1.49 


LOO 


2.95 


2.97 


.987 


6.40 


1.41 


2.44 


0.969 


2.74 


2.77 


2.94 


5.66 8 


3.18 


6.12 


.999 


2.81 


2.86 


4.83 


5.25 5 


5.13 


9.92 


1.01 


2.60 


2.66 


9.78 


4.73 8 . 


9.95 


19.6 


1.01 


2.53 


2.63 


10.1 


4.71 8 


10.3 


20.4 


1.00 


2.56 


2.66 


19.7 


4.27 5 


19.6 


39.3 


0.999 


2.75 


2.90 


49.5 


3.72 


47.2 


96.8 


0.997 


3.33 


3.46 



Table 3. Solubility of DCPD at 15 ° C 

Unadjusted quantities 





Compos 


ition of eqi 


iil. solns. 




k; xio 7 


Initial 








Ca 






*cid, Po 








P-Po 






Mx 10 3 


pH 


Ca 
Mx 10 3 


P 

Mx 10 8 




Davies 


D.-H. 


0.0967 


7.62 7 


0.736 


0.842 


0.988 


2.93 


2.95 


.499 


6.79 :J 


1.01 


1.51 


.998 


2.83 


2.86 


1.01 


6.36« 


1.41 


2.41 


1.01 


2.78 


2.81 


2.95 


5.60« 


3.20 


6.03 


1.04 


2.69 


2.74 


4.83 


5.20 9 


5.11 


9.96 


0.996 


2.59 


2.66 


5.04 


5.17 5 


5.19 


10.1 


1.02 


2.47 


2.54 


10.0 


4.680 


9.97 


19.8 


1.02 


2.48 


2.58 


20.1 


4.215 


19.4 


39.0 


1.02 


2.61 


2.75 


50.0 


3.62 8 


46.7 


96.4 


1.01 


2.89 


3.01 



Table 4. Solubility of DCPD at 25 °C 

Unadjusted quantities 



' Initial 


Composition of equil. solns. 


Ca 
P-Po 


K[ DC x 10 7 


acid, Po 
MX 10 3 


pH 


Ca 
MX10 3 


P 

MX10 3 


Davies 


D.-H. 


0.101 
.503 
1.02 
2.98 
4.91 
9.70 

19.7 

49.5 


7.58 
6.77 5 
6.33 
5.58 3 
5.I84 
4.72 8 
4.22 1 
3.65 2 


0.760 

1.03 

1.41 

3.27 

5.19 

10.1 

19.7 

47.2 


0.847 
1.49 
2.42 
6.21 

10.2 

19.9 

39.1 

97.6 


1.02 
1.05 
1.00 
1.01 
0.978 

.985 
1.01 

.983 


2.98 
2.89 
2.75 
2.84 
2.70 
2.98 
2.84 
3.29 


3.00 
2.92 
2.78 
2.90 
2.78 
3.10 
3.00 
3.44 



values of the Ca : (P— Po) ratio, which do not depart 
significantly from unity. In columns 6 and 7 are given 
the apparent solubility products, K' DC , (assuming no 
ion pairs) of the salt calculated with the use of the 
activity coefficients of Davies and of Debye-Huckel, 
respectively. It is apparent that the results exhibit the 
trend noted in previous work [1]: a minimum value of 



465 



Table 5. Solubility of DCPD at $7.5 °C 

Unadjusted quantities 



Table 8. Solubility of DCPD at 25 °C 

Adjusted quantities 





Composition of equil. solns. 




K' DC x 10' 


Initial 






Ca 
P — Po 






acid, Po 












MX10 3 


pH 


Ca 
MxlO 3 


P 

M X 10 3 




Davies 


D.-H. 


0.498 


6.75 5 


1.02 


1.50 


1.02 


2.84 


2.86 


1.00 


6.30 5 


1.41 


2.40 


1.01 


2.66 


2.69 


1.00 


6.29 9 


1.37 


2.43 


0.963 


2.60 


2.63 


5.02 


5.13 2 


5.29 


10.1 


1.04 


2.47 


2.55 


5.02 


5.12 


5.32 


10.2 


1.04 


2.43 


2.50 


9.99 


4.58 7 


10.5 


20.3 


1.01 


2.31 


2.42 


20.2 


4.13 7 


19.7 


40.3 


.980 


2.43 


2.58 


20.2 


4.13 


20.2 


40.4 


1.00 


2.43 


2.59 


49.8 


3.56 


48.9 


97.9 


1.02 


2.73 


2.87 



K' DC occurs near a pH of 5. The K^ c values in columns 
6 and 7 tend to converge in the more dilute solutions, 
since both expressions for ionic activity coefficients 
approach the limiting law in this region. 

Tables 6 through 9 summarize the data adjustments 
for 5, 15, 25, and 37.5 °C, respectively, as obtained 
under the conditions of Case III— C in section 3.5b 

Table 6. Solubility of DCPD at 5 °C 

Adjusted quantities 





Composition of equil. solns. 








Initial 








Ca 






acid, P 
Mx 10 3 










pH 


Ca 


P 


P-Po 


jLLX 10 3 


K DC X10 7 






Mx 10 3 


M x 10 3 








0.0984 


7.66, 


0.753 


0.851 


1.00 


2.49 


2.33 8 


.482 


6.84 4 


1.00 


1.49 


1.00 


3.10 


2.33., 


.996 


6.39, 


1.41 


2.40 


1.00 


4.18 


2.33., 


2.91 


5.62- 


3.18 


6.10 


1.00 


9.17 


2.34 8 


4.85 


5.25, 


5.07 


9.92 


1.00 


14.5 


2.33, 


9.78 


4.75, 


9.92 


19.7 


1.00 


27.7 


2.33., 


10.1 


4.72« 


10.3 


20.4 


1.00 


28.6 


2.33., 


19.7 


4.27« 


19.6 


39.3 


1.00 


52.8 


2.33 8 


49.2 


3.70., 


47.6 


96.8 


1.00 


119. 


2.34 



Table 7. Solubility of DCPD at 15 °C 

Adjusted quantities 





Composition of equil. solns. 








Initial 








Ca 
P-Po 


jJL x 10 3 




acid, P 








K DC x 10 7 


Mx 10 3 


pH 


Ca 
M X 10 3 


P 

M x 10 3 








0.0969 


7.63 6 


0.746 


0.843 


1.00 


2.61 


2.53 5 


.502 


6.79 5 


1.01 


1.51 


1.00 


3.25 


2.53.-, 


1.01 


6.35, 


1.41 


2.42 


1.00 


4.32 


2.53-, 


2.93 


5.60,, 


3.17 


6.10 


1.00 


9.33 


2.53-, 


4.88 


5.22 4 


5.07 


9.94 


1.00 


14.7 


2.53 5 


5.03 


5.20, 


5.21 


10.2 


1.00 


15.2 


2.54, 


9.98 


4.70, 


10.1 


20.0 


1.00 


28.8 


2.54, 


19.8 


4.23« 


19.7 


39.5 


1.00 


54.8 


2.54 


49.6 


3.65 


47.5 


97.1 


1.00 


126. 


2.54, 





Composition of equil. solns. 








Initial 








Ca 
P-Po 


/LtXlO 3 




acid, Po 
MX10 3 


pH 


Ca 


P 


K DC X10 7 






A/X10 3 


Mx 10 3 








0.101 


7.58 4 


0.751 


0.852 


1.00 


2.64 


2.56, 


.499 


6.76 8 


1.01 


1.51 


1.00 


3.27 


2.56, 


1.01 


6.32 5 


1.41 


2.43 


1.00 


4.36 


2.56, 


2.98 


5.55 8 


3.22 


6.20 


1.00 


9.52 


2.56 7 


4.97 


5.17 9 


5.15 


10.1 


1.00 


15.1 


2.56 2 


9.79 


4.69 


9.88 


19.7 


LOO 


28.6 


2.57 8 


19.7 


4.20 4 


19.5 


39.2 


1.00 


55.2 


2.56 4 


49.5 


3.61, 


47.1 


96.6 


1.00 


128. 


2.57 5 



Table 9. 



Solubility of DCPD at 37.5 °C 

Adjusted quantities 





Composition of equil. solns. 








Initial 






Ca 
P-Po 


/jlX 10 3 




acid, P 
MX10 3 


pH 


Ca 


P 


K DC x 10 7 






MX10 3 


Mx 10 3 








0.492 


6.73 7 


1.01 


1.51 


1.00 


3.18 


2.37 5 


.995 


6.29 7 


1.41 


2.40 


1.00 


4.25 


2.37 5 


1.00 


6.29 2 


1.42 


2.42 


1.00 


4.26 


2.37 4 


5.00 


5.13, 


5.20 


10.2 


1.00 


15.2 


2.37 4 


5.02 


5.12 8 


5.22 


10.2 


1.00 


15.2 


2.37 4 


10.1 


4.62, 


10.2 


20.3 


1.00 


29.6 


2.37 8 


20.3 


4.14, 


20.0 


40.3 


1.00 


56.9 


2.37 5 


20.3 


4.14 


20.0 


40.4 


1.00 


57.1 


2.37 5 


50.3 


3.56, 


47.2 


97.5 


1.00 


130. 


2.37 4 



using D.-H. activity coefficients. Comparison with 
tables 2 through 5 shows that in all systems the magni- 
tudes of the individual residuals are small relative to 
the values of the corresponding observables. The 
small size of the relative adjustments, (residual/ob- 
served value) X 100, is crucial, since the linearization 
of the condition functions, neglecting higher order 
terms, can be justified only if the above ratios are 
indeed small. 

The relative dispersion of the adjusted values of 
Ca and P, cr{c7f}/[/f, where the dispersions are given 
by eq (26), are not shown in the tables, but were well 
below the assigned relative errors given in section 2.3; 
likewise, the dispersion of the adjusted pH values were 
less than the absolute error given in section 2.4. Thus, 
the calculated dispersions support the reasonableness 
of the weights assigned in section 3.4. 

The last two columns of tables 6 to 9 give the ionic 
strength and solubility product, K DC , calculated with 
the adjusted data and final values of the stability 
constants. The values of K DC are given with more 
figures than are significant to indicate the adequacy 
of the adjustment procedure. Since K DC is defined 
(sec. 3.2) as the function (Ca — r)(P — r)lg Ca N, its near 
invariance over the range of concentrations supports 
the correctness of the model involving two ion pairs. 
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3.2 3.6 4.0 4.4 



5.2 5.6 6.0 6.4 6-8 "7.2 7.6 8.0 



FIGURE 1. Calculated solubility isotherms of DCPD at 5, 25, and 
37.5 °C. 

1 'ointS represent experimental data. 



In this respect it may be observed that the. variations 
in the listed K lK < values are clearly independent of pH, 
in contrast with those for K[ )C in tables 2 to 5. Finally, 
it may be observed that the K DC values are practically 
identical to the value of the solubility parameter K s 
listed in table 12. This is further evidence as to the 
soundness of the adjustment and the applicability of 
the model. 



4.2. Solubility Isotherms 

In figure 1 are shown the solubility isotherms for 
DCPD in terms of total concentration of calcium as a 
function of pH for three temperatures. These curves 
were computed, as explained in section 5.1, on the 
basis of the constants derived here. Thus, these curves 
are the loci of the adjusted data. Experimental points 
are shown on the same plot. It is apparent that the 
experimental data at 5 and 37.5 °C are well represented 
by the calculated isotherms. The reason for the less 
satisfactory fit for the data at 25 °C will be discussed 
in section 5.1. The isotherm for 15 °C lies very close 
to that for 25 °C and, for the sake of clarity, is not 
shown. The relative positions of the isotherms in 
figure 1, show that the solubility of DCPD decreases 
with increase in temperature. This behavior is similar 
to that displayed [14] by CaHP0 4 although the decrease 
in solubility is more marked in the latter case. Similarly 
smooth curves are obtained if the concentration of 
phosphorus is plotted against pH or if the concentra- 
tion of calcium is plotted against the concentration 
of phosphorus. 



4.3. Ion pairs 

The concentration of the ion pair [Caf^POj], F + , 
and the percent of bound calcium (100 r/Ca) for the 
adjusted data are given for all systems in tables 10 
and 11. In table 10 the results for 25 °C are shown for 
the calculations using the Davies and the D.-H. ionic 
activity coefficients. Only the results obtained with 
the latter coefficients are given in table 1 1 for the other 
three temperatures, since, as indicated by the values in 
table 10, the results using the two types of coefficients 
do not differ markedly. The concentration of the ion pair 
[CaHPOli], X l \ wa s nearly invariant in all systems at 
each temperature, its value being calculated in each 
cycle from eq (9), involving the quantities R and r, 
both of which vary with pH. This constancy was not 
assumed ab initio but is a consequence of saturation: 
eq (6) gives X°= K x ' K DC and therefore X° should be 
constant under the condition of saturation with 
respect to DCPD. The values of A" and its dispersion, 
given at the foot of each table, were computed after 
the final cycle, with K s replacing K DC in eq (6). 

4.4. Equilibrium Constants 

In table 12 are assembled the parameters A m u of 
eq (16) and the sets of values of K 8 , K x , and K y , and 
their respective errors at each of the four tempera- 
tures. (The standard deviation is discussed in section 
5a.) The values obtained using the D.-H. and the Davies 
activity coefficients are both shown. The differences in 
the values so obtained persists in a uniform manner 
from one temperature to another. Values o£K 8 , K x , and 
K y , given by eq (16), are plotted in figure 2; their stand- 
ard errors are indicated at the four experimental tem- 
peratures. The presence of a minimum in K x is in 
agreement with the behavior observed in many other 
weak electrolytes; usually they exhibit ionization 
maxima in the region to 60 °C and it has been sug- 
gested [15] that this is a property of all weak 
electrolytes. It is noteworthy that such a maximum 
(minimum for association) has been reported [16] 
for the second ionization constant of phosphoric 
acid. Furthermore, the maximum in K s can be antici- 

Table 10. Concentrations of ion pairs and bound calcium at 25 °C 



Debye-Hiickel 


Davies 


pH 


Y\MX10* 


100 r/Ca 


pH 


Y + , Mx UP 


100 r/Ca 


7.58 4 


0.0006 


8.76 


7.58, 


0.0004 


8.96 


6.76* 


.0037 


6.83 


6.76- 


.0025 


6.89 


6.32-, 


.0104 


5.35 


6.32 4 


.0069 


5.26 


5.55* 


.0631 


3.99 


5.55* 


.0405 


3.47 


5.17c 


.154 


4.25 


5.18, 


.0964 


3.41 


4.69,, 


.496 


5.68 


4.69 4 


.297 


4.21 


4.20 4 


1.59 


8.48 


4.20, 


.905 


6.05 


3.61, 


6.70 


14.4 


3.60* 


4.75 


10.2 


X° = 0.0652 ±0.0080 


X° = 


0.0669 ±0.00^ 


51 
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Table 11. 


Concentrations of ion pairs and bound calcium calculated with D 


-H. activity coefficients a 


5°C 


15 °C 


37.5 °C 


pH 


Y+,MX10 3 


100 r/Ca 


pH 


Y + ,Mx 10 3 


100 r/Ca 


P H 


Y+,Mxi0 3 


100 r/Ca 


7.66, 


0.0010 


14.0 


7.63 6 


0.0007 


9.34 


6.73 7 


0.0032 


9.72 


6.84 4 


.0069 


11.1 


6.79 5 


.0049 


7.31 


6.29 7 


.0088 


7.39 


6.39, 


.0197 


8.80 


6.35 9 


.0134 


5.85 


6.29 2 


.0089 


7.36 


5.62 7 


.118 


7.00 


5.60 


.0792 


4.67 


5.13! 


.137 


4.46 


5.25, 


.285 


7.68 


' 5.22 4 


.193 


5.16 


5.12 8 


.138 


4.46 


4.75, 


.937 


10.5 


5.20 3 


.203 


5.22 


4.62! 


.462 


5.45 


4.72 6 


.995 


10.7 


4.70 9 


.658 


7.22 


4.14 2 


1.46 


7.78 


4.27 6 


2.92 


15.4 


4.23 6 


2.05 


10.8 


4.14 


1.47 


7.79 


3.70.3 


11.7 


24.9 


3.65 


8.49 


18.0 


3.56 3 


5.97 


12.8 


X° = C 


.104±0.010 




X°= 0.0690 ±0.0079 


Z°= 0.0954 ±0.0153 



7 The values are not comparable across the table due to differences in initial H3PO4 concentrations. 

TABLE 12. Temperature dependence of equilibrium constants 



\nK m = 2 A mk T~ 2 + k [eq(16)] 



Debye-Hiickel 


Davies 




Amk 


Amk 




rn^^k 


1 


2 


3 


1 


2 


3 


1( 5 ) 


- 8403.5 
51090 
19373 


41.863 
-341.14 
- 122.81 


- 0.096779 
.58804 
.19943 


- 9418.4 
54735 
31465 


48.746 

- 365.90 

- 202.88 


-0.10848 


2{x) 


.63022 


3(y) 


.33085 






K m """""""" — ___ ___ ^C 


5 


15 


25 


37.5 


5 


15 


25 


37.5 


K,X 10 7 

K x 

K y 


2.33 8 
± 0.058 
446 
±54 
10.1 
± 1.7 


2.53 5 

± 0.044 
272 
±35 
6.62 
±1.23 


2.56, 
± 0.048 
255 
±35 

5.10 
±1.03 


2.37 5 

± 0.057 
401 

±71 
4.51 
± 1.44 


2.28 9 
± 0.058 
479 

±57 
9.69 
±1.83 


2.50o 
± 0.044 
283 
±36 

5.23 
± 1.29 


2.53 4 
± 0.048 
264 
±36 

3.67 
±0.98 


2.32 8 
± 0.057 
431 
±74 

3.28 
±1.48 



pated, as shown in section 5a, on the basis of the 
thermodynamic properties of DCPD itself. 

4.5. Thermodynamic Quantities 

The results of the thermodynamic computations are 
given in table 13. The values of the thermodynamic 
functions were calculated from eq (16) and the follow- 
ing equations: 

AG° = -RT\nK m = -R(A mi +A m J + A m ,T 2 ),(27) 

AS° = - (d&GoidT)=R(A m2 + 2A m . s T), (28) 
A// =AG + rAS° = /?(-^ mi + ^m,r 2 ), (29) 

where R is the gas constant. For DCPD (ra=l) the 
functions refer to the dissolution of 1 mole of CaHP04" 
2H 2 according to the reaction 

CaHP0 4 • 2H 2 = Ca++ + HPO= + 2H 2 0. (30) 



In the case of the ion pairs, [CaHPO°] (m = 2) and 
[CaH 2 P0 4 + ] (m = 3), each function applies to the for- 
mation of 1 mole of the ion pair from its constituent 
calcium and ortho-phosphate ions. The errors were 
computed by use of the appropriate moment matrix 
of the Amk, eq (23), and the derivatives of eqs (27-29). 
Whereas the errors for the quantities associated with 
the dissolution of DCPD, and formation of [CaHPO°] 
are within acceptable limits, this is not the case for 
those related to the formation of [CaH 2 POj]. 

The negative values of A S° for formation of the ion 
pairs at the lower temperatures are atypical. It would 
be desirable to check these values by another experi- 
mental method. On the other hand, the value of 
ASg s , -31.5 ±0.8 cal/mol.deg, for eq (30) is in 
satisfactory agreement with the value —32.56 calcu- 
lated from the S° (DCPD) [21] and the entropies [22] 2 
of the products in that equation. 



2 The revised value of — 12.7 cal/mol deg forS'is (Ca ++ ) was supplied by Dr. V. B. Parker, 
Thermochemistry Section, NBS. 
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2.8 




FIGURE 2. Temperature dependence (eq {16)) of the solubility product 
of DCPD and the stability constants for ion pairs. 

Length of vertical lines at the four experimental temperatures corresponds to twice the 
standard errors. 



It may be observed that the thermodynamic values 
appear to be relatively insensitive to the choice of 
activity-coefficient expression, and that the two sets 
of values differ by less than their errors. (The errors 
are approximately the same for the two sets.) By way 
of contrast, the large temperature dependence of the 
equilibrium constants produces marked variations in 
A//° and AS ; this is so for the values pertaining to 
the dissolution process even though the temperature 
coefficient of K s is relatively small. It is clear that a 
high degree of confidence in the results can be attained 
only by measurements made at short temperature 
intervals and extended over the widest practical 
temperature range. Thus, measurements on DCPD at 
other temperatures are highly desirable. 

5. Discussion 

5.1. Solubility isotherms 

The measurements reported in this paper were 
made on solutions that had reached equilibrium with 
respect to DCPD from undersaturation; the experi- 
mental design was selected on the basis of previous 
experience [1, 2] to assure attainment of equilibrium 
conditions. The consistency of the data, the smooth- 



Table 13. Thermodynamic quantities 





Dissolution of DCPD 


Formation of [CaHPO°J 


Formation of [CaH 2 PO+] 


t°C 


AC 

kj/mol 


A//° 

kj/mol 


J/mol.K 


AG 

kj/mol 


A//° 

kj/mol 


AS° 
J/mol.K 


AG 

kj/mol 


A// 

kj/mol 


AS° 
J/moLK 


5 

15 
25 
37.5 


«35.3±0.04 
6 35.3 

36.4 ±0.04 
36.4 

37.6 ±0.04 

37.7 

39.4 ±0.04 
39.4 


7.6 ±2.5 
8.5 

3.0± 1.2 
3.4 

-1.7 ±1.0 

-1.9 

-7.8±2.6 
-8.7 


-99.6 ±8.8 
-96.6 

-115±4 
-115 

- 132 ±3 

- 133 

-152±9 
-155 


-14.1 ±0.3 
-14.3 

-13.4±0.3 
-13.5 

-13.7 ±0.3 
-13.8 

-15.5±0.5 
-15.7 


-46±15 
-50 

-19±7 
-20 

10±7 
11 

46±19 
50 


-117±54 
-125 

-19 ±24 
-23 

79 ±23 
84 

201 ±63 
213 


-5.4 ±0.4 
-5.3 

-4.5 ±0.5 
-4.0 

-4.1 ±0.5 

-3.2 

-3.9±0.8 
-3.0 


-33±22 
-50 

-23±10 
-33 

-14±12 

-17 

-1±30 
4 


-100 ±79 
-155 

-67±34 
-100 

-32±40 
-46 

9±100 
23 



D.-H. activity coefficients. 



6 Davies activity coefficients. 



ness of the experimental plots, and the success of 
the adjustment procedure in the calculation of a solu- 
bility product constant applicable to all the systems 
(over a wide range of solution compositions) constitute 
further evidence that a true equilibrium was attained 
in the experiments reported here. The adjustment of 
the data yields points which He on the calculated 
isotherms at the several temperatures (fig. 1), and it 
may be seen from tables 2-9 that the only residuals 
that exceed 2 percent of the observed value are those 
of calcium, one at 25 and two at 37.5 °C. Furthermore, 
the residuals of positive sign constitute approximately 
half of the total at each temperature. 



The discrepancy at low pH values between the data 
at 25 °C and the isotherm computed for that tempera- 
ture can be attributed to the difference between two 
estimates of the value of K y . The adjustment of these 
data alone (Case I, sec. 3.5b) resulted in values for K s , 
K x , and K y of 2.59 X 10~ 7 , 244, and 8.2, respectively. 
Isotherms computed with these values, and with the 
values reported for K s and K x (table 12) but with 8.2 for 
Ky, both fitted the data in a satisfactory manner. Never- 
theless, we have elected to accept the lower estimate 
(5.1) of the value of K y for the sake of increased confi- 
dence in the values of K s and K x resulting from the 
adjustment of the combined data. 
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Thus, the curves in figure 1 constitute solubility 
isotherms for DCPD at the indicated temperatures; 
they represent solutions metastable with respect to 
CaHP0 4 [14] over all the experimental range, and with 
respect to hydroxyapatite [1], /3-tricalcium phosphate 
[17], and octacalcium phosphate [18] in parts of that 
range. Metastability, however, does not preclude the 
application of thermodynamic concepts to these sys- 
tems so long as precipitation of more insoluble calcium 
phosphates is impossible to detect during the time of 
the experiment. In this connection, if more basic cal- 
cium phosphates had precipitated, the apparent 
Ca:(P — P ) ratio for the dissolution reaction would be 
significantly less than unity; the ratios reported here do 
not depart significantly from unity. For example, at 
25 °C, table 4, the dissolution ratios may be weighted by 
propagation on the quantity Ca/(P — P ). The weighted 
average ratio obtained is 1.0067 with standard error of 
the mean 0.0087. With 7 degrees of freedom, the excess 
over unity of the average ratio would have to exceed 
0.020 to be significant at the 95 percent confidence 
level. Similar conclusions can be drawn from the data 
at other temperatures. Furthermore, with one excep- 
tion [1], petrographic examination of the solid phase 
after the saturated solution was obtained failed to 
show any solid other than DCPD. 

Effect of Temperature on the Solubility of DCPD. 
Examination of the curves in figure 1 leads to the tenta- 
tive conclusion that the solubility of DCPD varies 
inversely with temperature. However, since the 
observed solubility differences are very small, it was 
considered desirable to ascertain the significance of 
these differences by the following statistical procedure. 
The isotherms in figure 1 were generated, using the 
constants K x , K y , and K s as found by the adjustment 
procedure (D.-H. activity coefficients), by solving 
the electroneutrality eq (14) for calcium concentrations 
corresponding to a selected set of pH values. The 
value for total phosphorus used in this equation, as 
well as the derivatives of calcium with respect to the 
three equilibrium constants, were obtained from the 
saturation condition, eq (13). The variance of calcium 
was calculated at each pH by use of the derivatives 
and the covariance matrix, eq (24), of the constants 
as given by the adjustment. For two temperatures, 
T\ < T 2 , Ca(ri) was compared with Ca(r 2 ) at the 
95 percent level by direct computation of critical 
values [19]. 

The results indicated, with one exception, that 
throughout the larger part (4.0-6.4) of the experi- 
mental pH range, the solubility of DCPD decreased 
significantly with increase in temperature. The 
exception to this statement is that there was no signif- 
icant difference between the solubility at 15 and at 
25 °C. The resulting plateau can be related to the 
variation of the equilibrium constants with temperature 
shown in figure 2. Between 5 °C and a temperature near 
25 °C, the decrease in K x should tend to depress the 
solubility of DCPD, whereas from that temperature 
up to 37 °C the reverse effect would be expected. The 
temperature variation of K s clearly should have the 



opposite effects on the solubility over the range 5 to 
37 °C. Thus, in the region where the extrema in the 
curves occur, a decreasing rate of change of solu- 
bility with temperature should be anticipated. The 
results, therefore, indicate that the decrease in the 
relative concentration of ion pairs (tables 10-11) is a 
dominant factor in the decrease of the solubility of 
DCPD in the region 5 to about 25 °C. Above this 
temperature, the dominant factor appears to be the 
decrease in the value of K s . The decrease in K y over 
the whole temperature range would reinforce the 
effect of K x in the lower, and that of K s in the upper, 
region. 

It is noteworthy that the anhydrous salt [14] does 
not display any solubility plateau, which possibly is 
due to the monotonic decrease of K s with temperature 
for that compound. The variation of K s with tempera- 
ture for both compounds can be anticipated on the 
basis of their respective thermodynamic properties, 
as discussed in section 5.2a. 

5.2. Equilibrium Constants 

a. Dependence Upon Conditions and Hypotheses 

The solubility product constant, K s , obtained by 
the adjustment procedure is a least-squares estimate 
and is thus associated with a minimum sum of squares, 
i.e., eq (A-15). With a given set of data, however, the 
actual value of q 2 at the minimum, as well as the 
value of K s itself, depend upon the condition functions 
that are satisfied by the systems and upon the hypothe- 
sis of ion pairs. This is illustrated in table 14. Because 
the number of condition functions and parameters 
vary across the columns of the table, the results are 
best compared by reference to the sixth line, which 
gives the standard deviation of the adjustment, eq. (22). 

It is clear from the table that the use of three 
condition functions has definite advantages over the 
use of only one condition function, whether ion-pairs 
are considered or not. Furthermore, it is evident that 
the minimum value of the standard deviation is ob- 
tained by the use of three condition functions with the 
hypothesis of ion pairs present. Thus, the best set of 
values for K s , K x , and K y are those obtained with 
the use of three condition functions. The argument 
for the presence of ion pairs is strengthened by this 
comparison. In a previous publication [1] we reported 
a set of values for K s , K x , and K y at 37.5 °C which 
differ from those reported here. Although the present 
calculations are based on a larger sample of data, 
these differences are primarily due to the method of 
calculation. The criteria used in our earlier work 
[1] were valid but the procedure to establish an un- 
equivocal minimum in the standard variance was not 
fully adequate. Treatment of the previously published 
data [1] by the adjustment procedure of the present 
paper yielded results in agreement with the values 
reported in table 12. 

b. Temperature Dependence of K. s 

The appearance of a maximum in the K s of DCPD in 
the neighborhood of 25 °C, as shown in figure 2, de- 
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Table 14. Effect of the number of condition functions and the ion pair hypothesis 

and the equilibrium constants 


on the estimates of the standard deviation 


Function(s) a 


1 


1 +2 


1 + 2 + 3 


Case b 


11- A 


III-A 


1I-B 


I II— H 


II-C 


III-C 


Ion pairs 


Absent 


Present 


Absent 


Present 


Absent 


Present 


Degrees of freedom 


32 


26 


67 


61 


102 


96 






q 2 X 10 ,B , eq (A-15) 


139 


57.1 


182 


78.8 


322 


107 


5 X 10 8 , eq (22) 


2.08 ± 0.26 
2.85 ± .05 


1.48 ± 0.21 
2.59 ± .09 

237 ± 77 
4.8 ± 1.7 


1.65 ±0.14 
2.84 ± .04 


1.14 ± 0.10 
2.52 ± .07 
297 ± 60 
6.0 ± 1.4 


1.78 ±0.12 
2.85 ± .04 


1.06+ 0.08 


K, x 10 7 , 25°C 


2.56 + .05 


K x , 25 °C 


255 + 35 


K y , 25 °C 








5.1 + 1.0 













1 -Solubility product, eq (13); 2-Electroneutrality, eq (14); 3-congruent dissolution, eq (15). 
Cases as given in the tabulation of sec. 3.5.3. 



serves special consideration, particularly since such 
a maximum does not appear [14] in the solubility prod- 
uct of the closely related compound :} CaHP0 4 . It is 
shown below that this behavior is a possible conse- 
quence of the thermodynamic properties of the 
respective solids. 

The van Hoffs equation provides an expression for 
the temperature coefficient of K s . If the enthalpy of 
solution, A//J), for reaction (30) can be given explicitly 
as a function of temperature, then the temperature at 
which K s has a maximum, if any, may be obtained as a 
root of that equation. The enthalpy of solution is given 
by 

AH» = AH»+f AC p dT. (31) 

The change in the heat capacity can be obtained as a 
function of temperature with the aid of some reason- 
able approximations. Furthermore, since A//g can be 
calculated at 25 °C from the corresponding value of 
K s , the entropy of DCPD [21], and tabulated thermo- 
dynamic quantities [22] for the products of reaction 
(30), the constant AT/ft in eq (31) can be evaluated. In 
the present calculations using the value of K s at 25 °C 
given in table 12, a value for A//J) of —0.73 Kcal mol -1 
is obtained. 

The heat capacity for Ca ++ (aq) was roughly esti- 
mated to be — 6 cal mol 1 . 4 We assumed that the heat 
capacity for HPO= (aq) does not differ appreciably 
from the isoelectronic species SOf [23]. The addi- 
tional assumption was made that the heat capacities 
of water and the aqueous ions did not change with 
temperature over the range considered. The tem- 
perature dependence of the heat capacity, Cp, of 
crystalline DCPD (and that for CaHP0 4 ) was ob- 
tained by a least squares fitting of an equation of 
standard form to the data of Egan and Wakefield 
[21, 24] , 5 with the results (cal mol -1 ) 



3 It may be pointed out that the salt CaS() 4 • 2H 2 0, which is isostructural with DCPD, 
does show a maximum in its solubility near 15 °C [20], although it can not be inferred from 
this that a maximum exists in its solubility product. 

4 Values of the necessary partial molar heat capacities at infinite dilution (cal mol -1 ) are 
(1) for CI", -32.6 [22]; (2) for CaCl 2 , -71, kindly estimated by Dr. Vivian B. Parker of the 
National Bureau of Standards, Thermochemistry Section, by a rough extrapolation from the 
available data. A check on the value of —6 for Ca ++ derived in this way was obtained by 
using the value —74.4 for BaCU given in table (8—4^3) of reference [26]; a value of — 9 
was estimated in this way. 

5 Only "normal" points in reference [24] included. 



C>(DCPD) = 7.582 + 0.1361T-2.217X10 4 7- 2 , 

(50-310 K); 

Cp(DCPA)= 10.590 + 0.055797 1 - 4.896 X 10 4 r~ 2 , 

(80-310 K). 

In eq (31) for DCPD, ACp can thus be replaced by 
the expression 

ACp = -47.58-0.1367 1 +2.217Xl0 4 7 1 - 2 

and A//° is then found to be 19.58 Kcal mol 1 . Then 
the van HofFs equation takes the form 



/? — In K, (DCPD) 

o 1 



22.2 19.6 

- + - 



0.0476 



T 3 



T z 



-6.81 X10- 5 . 



The pertinent root of this equation corresponds to a 
maximum at 17 °C, which is in satisfactory agreement 
with the plot in figure 2. A variation of ± 10 percent in 
the assumed value for the CpiHPO^ ) and the use of 
— 9 cal mol -1 for Cp(Ca ++ ) did not affect this result 
significantly. 

Similar calculations for DCPA indicated a maximum 
in K s , far below the experimental range (i.e., ca— 22 °C); 
although this finding is not physically meaningful, it 
is consistent with the observed [14] monotonic decrease 
in K s for this compound in the range 5 to 37 °C. 

In view of the approximations and assumptions used 
in the calculations, only the qualitative aspects of 
these results are to be emphasized. They show that 
the presence of a maximum in the K s of DCPD, which 
contrasts with the monotonic decrease in that of 
DCPA, is consistent with the results of calculations 
based on the known thermodynamic properties of 
these compounds. 

c. Association Constants and Bjerrum's Theory 

The relationship of the association constant and the 
distance, a, between the constituent ions in an ion 
pair was derived by Bjerrum [25] on the basis of 
Coulombic forces 
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K a 



with 



= 4*irN nz 1 z 2 \e 2 
~1000|_ DkT 



\z1Z2\e' 
aDkT 



Q(b), 



(32) 



Values for the function Q{b) have been calculated by 
various investigators and summarized by Harned and 
Owen [26]. The use in eq (32) of the extreme values of 
the 95 percent confidence interval for K x (25 °C) and 
K y {25 °C) reported here, yielded the ranges 3.8-4.5 A 
and 5.9-6.1 A for the a parameters of [CaHPO$] and 
[CaH 2 PO|], respectively. Larger values of a are ex- 
pected for [CaH 2 PO J] than for [CaHPO°] on the basis 
of the valences of the constituent ions in the two ion 
pairs. 

Considering the interatomic distances usually found 
in calcium phosphates [27] and assuming sharing of 
the P0 4 -tetrahedron edge by the calcium, the distance 
between the phosphorous and calcium atoms would be 
in the order of 3 A. The apparent consistency between 
this result and the a values calculated above should be 
taken cautiously because if the two ions were indeed 
that close, polarization effects, not considered in equa- 
tion (32), would undoubtedly play a role. It is pertinent 
to point out that in the investigation by Dunsmore and 
James [28] on the dissociation of [MgSOJ] (isoelec- 
tronic with [CaHPO^]), the results of their conductivity 
measurements could be explained on the basis of eq 
(32), and were consistent with a value of 5.1 A for a. 

5.3. Comparison With Reported Data 

Greenwald et al. [29] studied the formation of the 
ion pair [CaHPO^] in mixtures containing H3PO4, 
NaOH, and CaCl 2 . From their pH measurements with 
a glass electrode (presumably at 25 °C) and also with 
a hydrogen electrode at 25 °C they calculated insta- 
bility constants at various ionic strengths; in their 
calculations, ionic concentrations were used instead 
of activities. The equation used by Greenwald et al. 
[29], to describe the dependence of their values on 
ionic strength yields a stability constant K x = ?>\6 for 
infinite dilution, in reasonable agreement with the 
value K x (25 °C) =255 reported here. 

Davies and Hoyle [30] reported instability constants 
at 25 °C for systems containing sodium and potassium 
phosphate buffers saturated with respect to CaI0 3 ; 
their pH measurements were made colorimetricaliy. 
The range of values obtained by Davies and Hoyle, 
^(25°C) = 417 to 625 and ^(25°C) = 11.5 to 12.5, 
agree in order of magnitude with those reported here. 

Greenwald [31] used the measurements reported by 
Bjerrum [32] in the dilute range of the homogeneous 
system Ca(OH) 2 — H3PO4— H 2 at 18 °C to calculate 
values for the instability constant of [CaHPOJJ]; 
these varied from 1.01 to 9.55 X 10~ 3 . The average 
value, 2.77 X 10~ 3 , yields a K x {\% °C)= 1/(2.77 X 10" 3 ) 
= 361, agreeing in order of magnitude with the value 
reported here at 15 °C. Greenwald [31] also recalcu- 
lated the value of K x from the data of Davies and 



Hoyle [30]; for this purpose, he calculated the pH 
from the experimental buffer ratios; his values for K x 
thus obtained covered the range from 56 to 100. As 
in previous reports by the same author [29, 33], ionic 
concentrations instead of activities were used in these 
calculations. Therefore, only order of magnitude 
agreement with the values reported here is significant. 

Similar rough agreement is found with the result of 
Gosselin and Coghlan [34], K X (25°C) = 12, who 
also used ionic concentrations instead of activities 
in their calculations. 

Recently, Chughtai et al. [35] have reported stability 
constants for the species [CaH 2 PO+], [CaHPO°], 
and [CaPOj] calculated from pH measurements in 
the system KOH— NaOH— H3PO4— CaCl 2 — H 2 0. By 
working in a low pH region (where the concentrations 
of the neutral and negatively charged species are 
negligible) those investigators calculated K y (25°C) 
= 25.6; this value was then used in systems at higher 
pH values to calculate ^(25°C)= 548; reported 
values at 37 °C were £,,(37 °C) =31.9 and K X (37°C) 
= 681. 

In the calculations of Chughtai et al. [35], the equa- 
tion describing the condition of electrical neutrality 
was solved algebraically with the measured concentra- 
tions and pH values considered exempt of error. 
Calculations based on perfect compliance of the data 
to the electroneutrality condition are extremely sensi- 
tive to small experimental errors; this was an important 
factor in the adoption of the statistical procedure 
reported in this paper. Furthermore, by the procedure 
of Chughtai et al. [35], the two constants are obtained 
as ostensibly independent quantities, whereas in 
fact, a strong correlation should be expected. 

The use in eq (32) of the values for K x and K y at 
25 °C reported by Chughtai et al. yields values of 3.1 A 
and 2.3A for the interionic distance parameters for 
[CaHPOJJ] and [CaH 2 PO + ], respectively. This is an 
unreasonable result in view of the crystallographic 
radii of the ions involved [27], and, furthermore, con- 
tradicts the expectation a[CaH 2 PO+] > a[CaHPO°] 
(sec. 5.2c). 

6. Appendix 

6.1. Generalized Least Squares for Independent 
Systems 

In this paper attention is focused on individual 
systems and this is reflected in the notation of table 1 
by the use of the superscript a to denote the ath 
system. Apart from this feature, our notation differs 
from that of reference [4] in the following respects: 
(1) n and r designate here the number of components 
in the ath observation and condition function vector, 
respectively, with a=l, . . ., N; in table [11.3.3], 
n and r refer to the total of such components for N 
systems. (2) The correction vector 6 with components 
#„ replaces — b. (3) The vector of adjustable constants 
(parameters) is here denoted by ft with initial values j8°. 
The definitions of reference [4] are inconsistent. The 
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equation 6j=6?-h&j, following [11.3, 1], uses bj and 
bj like fi v and /3£, with — bj^-d v . Elsewhere 6° is con- 
sistently used like /3°, but b appears in table [11.3.3] 
and subsequent equations as the correction vector. 
(4) The matrices J a and F a are the transposes of 
submatrices of J and F. Transposed matrices are 
indicated by the tilde. 



a. Solution of the General Normal Equations 

In order to determine the specific covariance 
matrices needed in these calculations it is necessary 
to identify the submatrices of x^ [11-3.3, 2]. The 
system of eqs [11.3.3.1, la] adapted to the analysis 
used in this paper can be exhibited in diagrammatic 
form as follows: 



V 









X 1 



I 



f 



L" 



k a 



f a 



J 



J a 



L N 


J" 


J N 






The matrix 4> is readily inverted by the method of 
submatrices [36] with the result 



1 [x pa ] -c- 1 ' 



where 



(A-l) 



i.e.,cDx|^=, 



b. Covariance Matrices 



P 



o 



Observations. The observations Yf are independent 
for all a and i, and therefore 



&(8Y°8Y?)=s*w?-% h 



and 



x «e = L" l bae-L"- x J«C-*pLr 



(A-2) 



Gov (y«, y$) = g ? (8y a 8y^) = s 2 W a l 8<#. (A-6) 



X p<* = C- 1 J a L a , 



(A-3) 



Conditions functions. The analogue of [11.3.3.1, 5b] for 

and C is defined by eq (18). In eq (A-2) 8 a(3 is the a P air of vectors is 

Kronecker delta. Since the submatrices of x rr an d 

x pr are x ap an< j x pa^ reS p e ctively, solutions analogous Cov (/°, f) = s 2 F a W a ~ 1 FP8 a p = s 2 L a 8 a p. (A-l) 

to [11.3.3.1, lc] may be written as linear combinations 

of the multivariate vectors/": Other covariance matrices. The following are sub- 

matrices of the indicated covariance matrices 

\a = ^ x a(if^ (A _ 4) d eriye( J in [4]. 

(1) Lagrangians [11.3.4, 2a] 
and 

Cov(\", \l i )=s 2 x ati . (A-8) 

e=2 $X P*fe. (A-5) 

The last equation, by virtue of the definition (A-3), (2) Residuals [11.3.4.1,2a] 

is just the solution of the reduced normal equations, 

eq(17). Cov(P\ V*)= #W *~ l FlC*F*W*~\ (A-9) 
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(3) Adjusted observations [11.3.4.2, 3b] 

Co\(U a , UP) =sW a ~ 1 8 ae - Cov (V, F' 3 )=5 2 w a ' s . 

(A- 10) 

c. Residuals and Sum of Squares 

The individual residuals vectors are subvectors of 
[11.3.4.1, 1]: 



and the weighted sum of squares is 

<? 2 = 2 V«W a V a . 



(A-ll) 



(A-12) 



This expression can be converted into a form more 
convenient for testing for convergence as follows: 
With the substitution of eq (A-ll) for V a and the 

subsequent identification L a = F a W a F a by eq (A-7), 
(A-12) becomes 



Now, with the aid of (A-2), (A-3) and (A-5), (A-4) 
for k a can be converted into 

k a = L«- 1 (f«-J«0). (A-14) 

With the substitution of (A-14) into (A-13) we get 
Q 2 = y Cf a -eh)L«-\f«-J«S); 

a 

carrying out the multiplication and using 

a 

eq (17), the final result is 

<? 2 = 2 ~f°L*- x f°-ece. (A-15) 

a 

6.2 Constants Used in the Calculations 



q 2 = ^ ^ a L a k a . 



The values for the temperature-dependent con- 
(A-13) stants are given in table Al. 



Table Al 



t°C 


Phosphoric Acid 


Water 
10 14 £ u ,[39] 


A(T)[4Q] 


B(T) [40] 


10 2 A,[37] 


10 7 A- 2 [16] 


10 12 k% 


5 


0.845 


0.523 


0.226 


0.1846 


0.4952 


0.3256 


15 


.782 


.587 


.320 


0.4505 


.5028 


.3273 


25 


.711 


.636 


.452 


1.008 


.5115 


.3291 


37.5 


.618 


.659 


.692 


2.485 


.5242 


.3318 



3 From the equation pk 3 = 12.45 - (t ° C - 18) X 0.015 of [38] . 



The values of the inter-atomic distance parameter 
a x used here are listed in table A2. 

Table A2 





H 2 P0 4 - 


HP0 4 = 


PO= 


H + 


OH- 


Ca+ + 


10 8 a,[41] 


4 


4 


4 


9 


3.5 


6 
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